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w I ND- TUN HEX INVEST I GAT I ON S 
ON ELEXURAL- TORSIONAL WING FLUTTER* 
By H. Voigt 

For the purpose of testing the theory of an oscillat- 
ing airfoil of two decrees of freedom, a wing was mounted 
in the wind tunnel "between two walls in such a manner that 
it could execute vertical (flexural) oscillations as well 
as torsional oscillations about an arbitrary axis parallel 
to the span. It was possible to vary the inertia and elas- 
ticity parameters and also to increase artificially the 
negligibly small natural damping of the system. The oscil- 
lations were recorded to a strongly magnified scale. The 
experimentally determined critical (or flutter) velocities 
fully a t :ree with the theoretical ones of Wagner and Kiissner 
within the limits of computational and measuring accuracy. 
An extremely narrow wing without end walls (three-dimen- 
sional problem) showed the same oscillations as one with 
end walls (two-dimensional problem). 



I , OBJECT OE THE TESTS 



In Germany in the past decade, methods have been pre- 
sented by Blenk and Liebers (reference 1) , by Kiissner (ref- 
erences 2 and 8), and by Wagner and Kassner (references 5 
to 7), for computing the critical velocity at which unstable 
wing oscillations may be expected (flutter velocity). Out- 
side of Germany the work of Theodorsen (reference 3) de- 
serves mention. Blenk and Liebers neglect the reaction of 
the vortices leaving the oscillating airfoil and hence give 
only useful approximations for oscillations with low re- 
duced frequency uu or high reduced velocity V. The ex- 
pressions of Kiissner, ffagner and Kassner, and Theodorsen, 
although derived from various initial assumption^, lead to 
the same solution determinants and give the same results. 
In the present paper the theoretical computat i on s are car- 

* "Unt er sue hung en von angef acht en Dr eh-Bi ege- Tragf lugel- 
schwingungen im ffindkanal . 11 Luf tf ahrtf or schung, vol. 14, 
no. 9," Sept. 20, 1937, pp. 427-433. 
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ried out partly according to Kussner, "but more largely ac- 
cording to the graphical method of Wagner and Kassner , 

The present tests are intended to verify the correct- 
ness of the tneory for the two degrees of freedom of wing 
deflection and wing rotation. For this piirpose, it was 
necessary to mount as light a wing as possible in such a 
manner that the natural damping of the system was a mini- 
mum, and so as to permit the parameters to he quickly and 
reliably varied. By far the greatest time expenditure was 
required for tue accurate determination of the system param- 
eters. By recording trie oscillations, an analysis of the 
phase and amplitude relations and an exact determination 
of the frequency were made possible. 

Tne tests extended over a period of 2-1/2 years. The 
results were mostly communicated to interested parties 
(ilenschel Airplane Company, Technical Flight Institute at 
Berlin, Institute for Aer o dy n&mis 3 at G-ott ingen , Junkers, 
and Heinkel) and were conducted in close cooperation with 
the first two named for -"hose encouragement and aid appre- 
ciation is hereby expressed. 

The present report includes only comparison with the 
theory within the scope mentioned. Results of first tests 
with ailerons, the relations between steady and unsteady 
oscillations (flutter), the effect of externally induced 
vibrations, and the effect of vibration dampers are dealt 
with in reference 4. 

II. COLLECTION OF SYMBOLS USED 



c, flexural stiffness of wing 
e t, position of elastic axis behind center of pressure 
T) t, elastic radius 



- 4. c 



i t, radius of gyration of wing with respect to center 
of gravity 

m *a 

• (jL , reduced wing mass or mass ratio = — - 

a L 

my, mass of wing with end a t tachmen t par t s , oscillat- 
ing springs, additional weights, and recording 
apparatus 
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mass of the air cylinder described by the wing 
chord 

torsional stiffness of the wing 
frequency of the oscillations 

position of the center of gravity "behind the 
elastic azi s 

wing chord 

wind velocity 

critical or flutter velocity 
reduced velocity 
reduced f r e quency 

phase "between flexural and torsional oscillations 
III. DESCHIPTION 0? THE TEST SET-UP 



The results, presented in sections VI and VII, were 
obtained on a wing of 250 mm chord, 800 mm span, and 30 mm 
thickness. The Wing had an HAGA profile 0013/63, was con- 
structed of wood, its forward portion covered by plywood 
and after portion by Japan tissue, and could he considered 
a s s t i f f • 

The wing was suspended between the two end walls, the 
suspension consisting of two end parts attached to the tips 
of the wing and to each of which was fixed four helical 
springs adjusted to the same stiffness. To obtain the drag 
a drag bracing member was employed, connected at each end 
to a pin at 32 percent chord. The position of the elastic 
axis and the elastic radius was determined by the points of 
attachment of the vertical springs and their distances 
apart. For varying the mass, the position of the center of 
gravity and the radius of gyration weights were employed, 
which were clamped to the end attachment parts (figs. 1, 2, 
and 3) ♦ 

The oscillations were recorded by two concave mirrors, 
one of which recorded the rotation (or torsion) of the wing 



m L , 

Ivi/cp, 

a t , 
t , 

v , 
V, 
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and the second, by means of a simple lever transmission, 
the deflection, i.e., tne vertical motion, the record be- 
ing made on light- sen si tive paper (fig. 2). This method 
has the advantage that the rotation need not necessarily 
be measured from two deflection amplitudes but can be re- 
corded at the same position of the film with the deflec- 
tion of an arbitrary point, for example, the elastic axis. 
It is possible, furthermore, to c noose the recording scale 
arbitrarily large. The scaLa was determined by calibra- 
tion, 1 mm on the record corresponding to 0.157 mm of wind 
deflection and 0.1315° rotation. 

The natural damping of the set-up was negligibly small 
as is snown by the oscillation curve (fig. 4), It was pos- 
sible to increase the damping artificially by permitting 
the lower springs, which were situated in vessels (fig. 3), 
to operate in glycerine-water mixtures of various composi- 
tion. The largest investigated damping decrement £ = 0,4 
lies above all value of the material damping decrements of 
the wing that occur in practice. 



IV. DETERMINATION OF THE PARAMETERS 



The vibrating mass of the system consists of the mass 
of the. wing, the side attachment 3 , the weights and the sus- 
pension with dra.g bracing and mirror mounting. The first- 
named components were determined by weighing , the latter by 
frequency determinations in oscillation tests. The associ- 
ated oscillating air mass was determined by oscillation 
cornputat ions. 



The position of the center of gravity of the win.?; with 
the side-attachment parts was determined by weighing of the 
wing suspended on two V wires. The center of gravity of 
the total system was computed. 

The elastic data, that is, flexural stiffness, tor- 
sional stiffness, and elastic ^axis , may readily be deter- 
mined theoretically from the stiffness- of the individual 
springs and their- p o s i t i on s ,o f a 1 1 achraen t , The relations 
were affected, however, by- the drag br-ac ing , the finite, 
vertical distance apart of the spring- at tachment points at 
the sides, the finite length of the spring extension, and, 
to a slight extent, by the lever of tne deflection mirror; 
so that for the accurate determination of the elastic data,, 
stiffness measurements with various spring distances and 
elastic axes were required. 



N.A.C.A. Technical iJeniorandurn Ho. 877 



The accurate determination of the moments of inertia 
was difficult, since the particularly light wing was af- 
fected during its oscillation by the mass of air vibrating 
with it. The most accurate method appeared to "be by means 
of the oscillation tests of the wings in the test set-up, 
the wing "being fixed at its axis through the center of 
gravity so that it could perforin only torsional vibrations, 
and the paper covering cut out and attached rolled up. 

The damping was likewise determined by oscillation 
tests. In order to exclude the quite considerable air 
damping, whose decrement $ /tt was determined up to 0.03, 
the wing was replaced by a tube which connected the two 
side attachments. The system was periodically excited by 
rubber tension. The rubber tension was removed by burn- 
ing a thread when the system in its flexural or torsional 
oscillations vibrated in resonance. In this way, natural 
vibrations without secondary vibrations were obtained. 
This method was also applied in the determination of the 
moments of inertia. 



Y. TEST PROCEDURE, 
ANALYSIS OF SEVERAL 0? IHB OSCILLATION RECORDS 



The wind velocity was first gradually increased in 
the wind tunnel until flutter was set up. This was gener- 
ally uniform so that the oscillation was of constant am- 
plitude as shown in figure 5a. To each wind velocity 
there corresponded a definite amplitude (fig. 6). In con- 
trast to the theory, which is based on the assumption of 
infinitely small amplitudes, a damping by the air thus oc- 
curred at finite amplitudes. The sharp theoretical sepa- 
ration line between the unstable and damped oscillations 
becomes a broad band of the stable region and the upper 
limit of the critical velocity is determined by the 
strength of the airplane wing. In the diagrams, figures 
3 to 13, velocities at which uniform amplitudes of about 
two percent of the chord occur are denoted by the sign 

The wind velocity was then slowly reduced and the 
point (>) noted at which the oscillations stopped of them- 
selves. This point must correspond most accurately to the 
theory, since here the limiting case of infinitely small 
oscillations with vortex field extending very far behind 
was most nearly satisfied. 
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In particular cases, especially for small distance of 
the center of gravity "behind the elastic axis and at high 
dynamic pressures, it was not possible to obtain stable 
oscillations; oscillations were more often set up suddenly, 
the wing striking against the stops and being brought to 
rest again. These cases, of which a mild example is plot- 
ted in figure 5b are denoted in the diagrams by the sign 
< >, The irregularity is explained by the fact that the 
air-stream disturbances originating from the blower, which 
could not be entirely removed by the honeycomb at some par- 
ticular instants, predominated over the action of the air 
forces. In still more extreme cases, the oscillations on 
account of their irregularity could hardly be considered as 
flutter oscillations (denoted by u) . 

Figure 5c shows, finally, a wing at the limit of stat- 
ic stability at which limit small changes in the angle of 
attack lead to an upward or downward tipping of the wing. 

Figure 7 shows oscillation records with the test ap- 
paratus in its first form in which by dry friction damping, 
whose decrement decreases with increasing amplitude, an un- 
stabilizing value was obtained, resulting in a particularly 
unstable oscillation. Frequencies, amplitudes, and phase 
angles between rotation and bending are given' oh figures 5a 
and 7. Further analysis of the curves and a detailed dis- 
cussion of the effect of the phase angle and the effect of 
external excitation on the phase angle are found in refer- 
ence 4. . 



VI. COMPARISON WITH THE THEORY FOR THE CASE WITHOUT DAMPING 



In order to test the correctness of the theory as com- 
prehensively as possible not only individual points but also 
entire test series were compared with the theory. In figure 
8 the position of the center of gravity, and hence also the 
radius of gyration, was varied . This results in only a slight 
change in the critical or flutter velocity. At the center 
of the curves there is a slight scattering of the computed 
and test points. For the computation the graphical method of 
Kassner was employed. In the tests a wire of the spring sus- 
pension appears to have ruboed somewhere. The test results 
are extremely sensitive to even the smallest disturbances of 
this kind, * 



In figures 9, 10, and 11, the elastic radius is the 
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principal variable. Variation of the latter also producos 
a slight variation in the radius of gyration, The three 
series of tests differ in the various values of the radius 
of gyration, the added weights "being of equal size and with 
equal position of the center of gravity but with various 
distances from the center of gravity of the entire system. 
In these series of tests the critical velocity is strongly 
affected by a change in the elastic radius. In spite of 
the extraordinary sensitivity of the measurement s t it is 
possible to speak of sufficient agreement between theory 
and experiment. A fact to be noted is that in the measure- 
ments flutter vibrations could "be observed up to the neigh- 
borhood of zero velocity. 

In figure 12 the results of three series of tests are 
simultaneously presented for which only the stiffness of the 
springs was varied. This results also in a slight change 
in the radius of gyration and the oscillating mass. The in- 
dependent variable in this series of tests is the position 
of'the center of gravity. In the three tests, there were 
equal reduced velocities V, that is, equal oscillation 
conditions for different wind velocities. Since, in each 
case there is agreement "between tneory and experiment, it 
is thus shown that a change in the Reynolds Number has no 
effect on the results. At the small Reynolds Numbers em- 
ployed, the M angle-of -at tack sensitivity" d c a /d a of the 
stationary flow is 25 percent lower than in the normal 
flight range as was established, in agreement with similar 
measurements, by investigation of the wing itself. Figure 
12 shows, as also the general agreement between theory and 
experiment in the previous figures, that during the oscil- 
lation all boundary layer processes are to a large extent 
eliminated and that it is correct to use in the computation 
the theoretical angle of attack sensitivity (abou.t 2 rr) . 

In order to test also very light wings with the least 
mass ratios (\i = 2,6), a wing of double the chord and equal 
in thickness ratio was tested on the same test set-up. It 
was here not possible to obtain complete agreement between 
theory and experiment.. On account of the turbulence of the 
jet,' that is, the irregularity of "the flow coming from the 
blower or guide vanes and which was not entirely eliminated 
by the honeycomb, the disturbance of the oscillatin airfoil 
was greater the smaller its mass. In the case of very light 
wings, the critical velocity can thus be considerably in- 
creased by these di s trubance s • It is also shown that the 
jet of a wind tunnel has a large tendency to oscillate. 
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Through the use of a wins; chord of 500 mm, the jet becomes 
highly loaded. It is therefore conceivable that to the 
oscillations of the wing there are added jet oscillations 
of some kind which may displace the critical velocity con- 
siderably "below or above. It is therefore planned to car- 
ry out several supplementary tests in the larger wind tun- 
nel on light wings and to provide special honeycombs for 
quieting the jet. From these tests, it may be concluded 
that similarly in flight through strongly turbulent air 
f lexural- tor sional flutter oscillations of particularly 
light airfoils may be delayed. 



VII. COMFARISOV STITH THE THEORY 
IN THE CASS OF VARIABLE DAMPING 



On figure 13 the critical velocity has been plotted 
against the variable damping decrement. The test data are 
taken from figure 8. Theory and experiment show satisfac- 
tory agreement so that also for the damping the theoretical 
expression should be considered as correct. 

On figure 14, finally are plo.ttei several results 
that have been obtained with the first test set-up in 1935. 
The natural damping was here still rather large and corre- 
sponded approximately to the material damping of the wing. 
The damping was determined by vibration tests and includes 
the air damping of the oscillating airfoil. Here, too, 
there is good agreement between theory and experiment. 

VIII. TESTS FOE THE THREE - l) I MEN.SIOSFAL PROBLEM 



The starting point was the following consideration. 
On a wing element which possesses a very small span in 
comparison with its chord, there is first i nv e s t i ga t ed ■ t he 
two-dimensional oscillating condition by placing the wing 
between two end walls. If these end walls are now removed, 
the wing element is exposed to a more intense flow at its 
tips than a wing element as part of an airfoil with, normal 
aspect ratio. All the effects of the tip flow must there- 
fore in this model appear particularly emphasized.. Con- 
versely, an agreement of the critical velocity in this lim- 
iting case would indicate that the effect of the. tip flow 
of the wing does not show up in the wing flutter. 



K .A. 3. A. Technical Memorandum No; 8 77 



The 



model employed had the following data (see fig. 



15) : 



Chord 



500 mm 



Aspect ratio 



0.5 



Span 



300 mm 



Flexural elasticity 0.58 kg/cm 



The vring was suspended on four springs, which were 
situated outside the air stream, four simple V wires leading 
from the springs to the wing. In addition there was provid- 
ed a longitudinal member for taking up the drag. The am- 
plitudes at the forward plane of the spring were read on a 
gage. The side walls were entirely independent of the test 
set-up. In order to test the effect of the wing without 
walls it was made porsiblc to add to the wing end section a 
rounded piece whose weight and radius of gyration were 
otherwise replaced in the wing "by added wioghts on the wing. 
The elastic radius was varied by cutting down the springs, 
and the mass and the center of gravity as well as the radius 
of gyration by adding weights to the wing. 

In all cases, it was found that the removal of the 
walls did not result in an increase in the critical veloc- 
ity "but occasionally even in a small decrease (table I). 
from this result it is to "be concluded that no reduction 
factor need be applied for the tip flow of the three-dimen- 
sional wing. The lowering of the critical velocity on re- 
moval of the walls is to he explained by the fact that in 
the case of the short and light wing the boundary layer be- 
tween the wing and wall disturbs the process and thereby 
appreciably raises the critical velocity. The critical 
velocities obtained without the walls are thus more reliable 
than those obtained with the walls. A rounding off of the 
wing-tip section had no appreciable effect on the final re- 
sult . . 

On table I are also presented the theoretically comput- 
ed critical velocities. The deviations, which are at times 
quite considerable, between theory and experiment arc prob- 
ably to be ascribed to the causes discussed above of tur- 
bulence and vibration of the jet. Bearing in mind, however, 
that the wing is so sensitive to even the smallest distur- 
bances the agreement in the two- and three-dimensional cases 
is surprisingly good. 
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IX. FURTHER TESTS PLANNED IN TEST PROGRAM 



As a continuation of the above tests the following 
further investigations are planned: 

1. Short supplementary tests on a particularly light 

wing in the moderate- si zed wind tunnel. 

2. Investigation of a symmetrical tapered wing of as- 

pect ratio 1:6 in the large wind tunnel to 
check the assumptions on the three-dimensional 
problem and to investigate the coupling "between 
the torsional and flexural vibrations. 



3. Further investigations with ailerons (ref. 4) for 
checking the theory with three degrees of free- 
dom and particularly for explaining the effect 
of the aileron slot on the oscillation process, 
an effect which is probably quite considerable 
and which so far has not been taken into account 
by any theory. It is planned to have the models 
so constructed as to make it possible to test 
individual industrial designs by means of models 
in the wind tunnel within the shortest time pos- 
sible and thus effect a time-saving in making 
the extremely extensive computations required 
for the case of the wing with aileron (three de- 
grees of freedom) . 



X. SUMMARY 



It is shown that the critical velocity obtained exper- 
imentally on a model wing suspended between two walls agreed 
with the value theoretically computed. For a wing of ex- 
tremely low aspect ratio the removal of the end walls does 
not raise the critical velocity. It is thus permissible to 
apply the computations of the two-dimensional theory without 
a correction factor to the three-dimensional case. 



In the case of extremely light wings of large chord, 
no absolute* complete agreement between theory and experi- 
ment could be obtained on account of the mutual interfer- 
ence effect between tunnel and wing. Further tests in this 
direction in another tunnel and additional investigations 
on wings with ailerons are planned. 

Translation by S. Reiss, 

National Advisory Committee for Aeronautics. 
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TABLE I: Compression Between Two- and Three-Dimensional Wing a) Wing oetween End Walls 
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Figs. 1.2.3.4.15 





Figure 1.- Wing with aide 

attachement parts. 




Figure 2.- Side attachment 

with added weight 
and mirror mounting. 




Figure 3.- 



Suspension of wing 
between end walls. 



i4 



flexure Torsion 




Figure 15.- Tests for three dimensional 

problem. Photo shows wing 
employed without side walls and with 
rounded edge pieces. 



Figure 4.- Portion of oscillation 

( flexural) curve for 
de terming the natural damping of 
the system. Wing replaced by tube. 
Oscillating mass l/g.240 g corre- 
sponding toU s 4.8 JL 0.0012. 
V = 55 per sec. 
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Figs. 5,6 



10° 


• .05 


8 - 


-.04 


6 1 


-.03 


4 ■ 


■ .02 


2 


.01 





a „ Flexure 



Torsion 



c 






x /v /v Vvwv] 





Figure 5.- Examples 

of observe! 
oscilliation curves, 
a) undamped, uniform 

oscillation 
pha st angle 1 00 0 
( denoted by sign = 
in figs, 8 to 13). 
b) 3&clllat:ioir 
disturbed by 
fi irregularities in jet 
)f\ ( lenoted by ( <>). 
c) Flexural- torsional 

oscilliation near 
the limit of static 
stability. Slow 
travel of the wing 
upwards . 
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Figure 6.- Increase 

in oscil- 
lation amplitude with 
velocity. 
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Figs. 7,12.13.14 




Figure 7.- Osci Illation record with determination 

of frequencies, amplitudes and phase* 
Damping for the most part dry friction. 
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Figure 13.- Effect of damping 

in theory and 
experiment. Data according to 
fig- 8. 
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Figs. 8,9.10.11 
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Figure 8,- comparison between Figure 9.- Comparison between compu- 

computation and test. tation and test. Elastic 

Position of center of gravity radius varied. Radii of gyration 

varied (series 13). small (series 21/25/2). 
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Figure 10.- Radii of gyration greater 
than in fig. 9 (series 
21/25/4) , 



than in fig. 9 and 10 
(21/25/6). 



